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We consider the possibility that the dark energy is made up of two or more independent compo- 
nents, each having a different equation of state. We fit the model with supernova and gamma-ray 
burst (GRB) data from resent observations, and use the Markov Chain Monte Carlo (MCMC) tech- 
nique to estimate the allowed parameter regions. We also use various model selection criteria to 
compare the two component model with the ACDM, one component dark energy model with static 
or variable w (XCDM), and with other multi-component models. We find that the two component 
models can give reasonably good fit to the current data. For some data sets, and depending some- 
what on the model selection criteria, the two component model can give better fit to the data than 
XCDM with static w and XCDM with variable w parameterized by w — wo + w a z/(l + z). 



I. INTRODUCTION 

Observations have shown that more than 2/3 of the 
total cosmic density is made of an unknown "dark en- 
ergy" , which cause the expansion of the Universe to ac- 
celerate [H0. The cosmological constant, with equation 
of state w = — 1, is the simplest form of dark energy, 
and it fits many observations reasonably well. How- 
ever, from a fundamental physics point of view, it suf- 
fers a fine-tuning problem of 122 order of magnitude. 
Furthermore, it appears to be a great coincidence for 
its density to be just comparable to the matter den- 
sity today 3- Many alternative models of dark en- 
ergy, e.g. quintessence SJjJE 0, 0, i, 13 , phantomfidj , 
auintom[12j. K-essence 

UMM, Chaplygin gas[H[ll[I3, 

and so on, have been proposed, but at present none of 
these is clearly superior than the cosmological constant. 
Most of these models consists only a single component of 
dark energy: either it is phenomenologically described by 
an equation of state parameter w : p = wp (w can vary, 
and many different forms of parameterization have been 
introduced), or by a dynamical field, although in some 
cases two scalar fields are introduced. 

In the present paper, we consider dark energy models 
which consist of two independent components with equa- 
tion of state w\ and W2 respectively. This is motivated 
by the following considerisions. (1) If there is indeed un- 
known physical laws which give rise dark energy, it may 
well allow more than one form of it to be present. (2) 
Furthermore, if one describe dark energy by scalar fields, 
then they do not necessarily have one single equation of 
state, especially in the case with multiple fields, it is quite 
plausible that for each field the correspondin g equat ion of 
states are different. (3) Some observations [l8|, [l9| seem 
to favor models with w < — 1 which necessitate negative 
kinetic terms (phantoms), or with variable equation of 
state in which w vary accross -1 (quintoms), although in 
more recent observations [M HH these trends are weak- 
ened. One may wonder if such unusual behavior could 
be more naturally explained by considering two dark en- 
ergy components, with different equation of states, and 
dominate at different redshifts. For example, in Ref. [Hj] 



two scalar fields were invoked to reproduce the phantom 
behavior. 

We model the two independent dark energy compo- 
nents with equation of state parameter w\,w%, which 
can either be fixed, or variable with redshift. We use 
the Markov Chain Monte Carlo (MCMC) technique to 
find the probablity distribution of the models in parame- 
ter space, which allows us to see how the two component 
dark energy model fares with observation. The data we 
used are supernova type la (SN la) and GRB data, which 
is currently the most direct and powerful probe of dark 
energy. In this analysis we have not included other con- 
straints , esp. the cosmic microwave background (CMB), 
as we do not expect the fixed effective equation of state 
description used here is really good at very high redshifts. 

To compare our model with single component dark 
matter models and see if there is indeed any evidence for 
two component of dark energy, we calculate the best fit 
X 2 , as well as the AIC, BIC and Bayesian evidence(BE). 
The models being compared include the ACDM, XCDM 
(phcnomclogical model of single component dark energy 
described by an equation of state), and our two compo- 
nent model XiX 2 CDMl. 

We have used several different data sets to perform 
our fit, e.g. the Gold04 [3 and Gold06 [3 data set, 
the SNLS [H data set, and the ESSENCE data set[H|. 
There are some differeneces in the result. To avoid repe- 
tition, we presented the result of only one data set in 
detail, this is a combination of supernovae and GRB 
data. The supernovae data consists 182 high-quality SN 
la samples 1 selected from the Gold06 [3, SNLS [if and 
ESSENCE U data sets, which includes 30 HST super- 
novae, 47 SNLS supernovae from Gold06, 60 ESSENCE 
supernovae, and 45 nearby supernovae from WV07 [23| . 
Two different light-curve fitters, MLCS2k2 [Mill, and 
SALT |20| , are used in these data sets. As they are con- 
sistent with each other [l!| [2l|, all of the SN la data 
used here are fitted using MLCS2k2 algorithm to avoid 



1 Although the number of sample coincides with that of Gold06, 
this is in fact a different data set. 
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normalization [27j . The GRB data provides a good com- 
plement to the SN la data [H, [H M, HH [H, M, H 11 , 
as the GRB can be observed at much higher redshift, 
thus providing a better baseline for detecting variation 
in the dark energy equation of state. The GRB data set 
is constituted of 27 GRB samples in Ref . [2{| , which are 
generated with the E peak — E 7 correlation discovered by 
Ghirlanda et al.(2004) and is one of the tightest correla- 
tion for GRB. We have not used all the 69 GRB samples 
in that paper, because there is larger uncertainty in the 
other correlations. The method of calculating /^grb is 
described in Ref. [2{|, but we only use the E pea k — E 1 
correlation in the paper. We marginalize Hq at last as 
we did for the supernovae case. 

The structure of this paper is as follows: in §2 we 
present the formalism and methods, including a brief 
summary of the MCMC approach for computing prob- 
ablity distribution of the parameters, and brief summary 
of model selection criteria. In §3, the results of our fit 
are shown. In §4, we will compare a few models. Finally, 
we summarize our results in §5. 



II. METHODS 
A. Model 

We consider a model of two independent components 
of dark energy, with static equation of state w± and W2 
respectively. The expansion rate H(z) is given by 

H 2 (z)=H 2 n(z;0), (1) 

where, in the flat case: 

n(z;0) - a mo (i + zf + n Xl (i + z)^ 1+w ^ 

+(i - n mo - n xi )(i + z )3(i+«>2) (2) 

and in the non-flat case: 

n(z;0) = {i - n mo - n Xl - n X2 ){i + z) 2 + n rno (i + z) 3 
+n Xl (i + z) 3 (!+^) + n X2 (i + z) 3{1+w ^ (3) 

Here 9 is the cosmological parameter set consisting of 

(Q mo ,fl Xl , wi,w 2 ) in flat case ^ 



(f2 m „ , n xi , fl X2 ,wi,W2) in non-flat case 
The luminosity distance can be written as 

d L (z;e) = (l + z)\n k \-hmn{\n k \ L * xJ'J^L}, (5) 

where sinn(x) = sinh(a;), x, sin(a;) for open, flat and close 
geometries respectively. 



B. Analysis 

Given a cosmological model defined by n parameters 
= (0i, ■ ■ ■ ,9 n ), and a data set consist of N quantities 



d = (di, . . . , d/v) with Gaussian-distributed errors <7d = 
(ai, . . . , cr/v), the likelyhood function can be written as: 



C(d\0) = -j^e- 
where for un-correlated data points, 



(6) 
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If both SN la data and GRB data are used, we have 



2 _ 2 i 2 

X — Xsn sc i + Xgrb • 



For the SN la data, 



XSN 



(*) = £ 



(Hobs(Zi) - Vth(Zi)) 2 



(8) 



(9) 



where the theoretical value of distance modulus fith(zi) 
is given by 

Vth(zi) = 5log 10 d L (zi) +25 

= 51og 10 ^ L (z l )-51og 10 ^q + 42.38, (10) 



and 



Dl(z) = — x d L (z). 
c 



(11) 



The Hubble constant Hq — 100/iokms _1 Mpc _1 may be 
marginalized over as a nuisance parameters, as described 
in [36j , and we would obtain 



xln s JO)=A(0) - 



B 2 (9) 
C ' 



(12) 



where 



A(0) 
B{ff) 
C 



N 

E 



JV 

E 

i=l 
N 



\Uob s (zi) - 51og 10 (D L (z i )) 



^obs(zi) - 51og 10 (D L (zi)) 



E rr2 



(13) 



Alternatively, we can integrate over Hq, and be left with 
an additional constant term ln((7/27r) which would not 
affect our results given a data set [H, [33, HI]. 



C. MCMC 

Given an observational data set d, the posterior dis- 
tribution of a parameter set 9 is, according to the Bayes 
theorem, 



p(0|d) = 



C(d\0) P (9) 
JC(d\0)p(9)d9'- 



(14) 
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where p(0) is the prior probability distribution, and 
C(d\8) is the likelihood for obtaining the data set d given 
the parameter set 9 and is given by Eq. [6] To obtain the 
posterior probability distribution of the parameters, we 
employ the MCMC technique to generate random sam- 
ples in the parameter space. This method has several 
advantages over grid-based approach. Most importantly, 
the computational time cost increases linearly with the 
number of parameters, so even for a large number of 
parameters the estimate can be done within an accept- 
able computation time. Additionally, since the MCMC 
method generate samples from the full posterior distribu- 
tion, it gives far more information than the marginalized 
distributions [39l lioL lilj] . and dose not require a Gaussian 
distribution of the likelyhood[I3]. 

With the MCMC technique, a chain of sample points 
which is distributed in the parameter space according to 
p(0|d) is generated by Monte Carlo. The Metropolis- 
Hastings algorithm with uniform prior probability distri- 
bution is used to decide whether to accept a new point 
into the chain by an acceptance probability: 



7n+l 



l*n) 



. j P(0n+l|d) q(g n |fln+l) 

1111 \ P (0 n |d) q(e n+1 \e a ) 

. f £(d|fl n+1 )q(fl n |fl n+1 ) 
1111 \ C(d\9 a ) q(0„+i|0„) 



,1 (15) 



where q(# n +i|#n) is the proposal density of proposing a 
new point 9 n+ i given a current point 9 n in the chain. 
If a = 1, the new point 9 n+ \ is accepted; otherwise, 
the new point is accepted with probability a. The trials 
are repeated until a new point is accepted, and then we 
set 6 n = 0n+i- In our computation, we set a uniform 
Gaussian-distributed proposal density for every point 
which is independent of the position on the chain, so 
that q(0n+i|#n) and q(0 n |0n+i) are cancelled, we then 
have 



7 n+l| 



£(d|0 n+1 ) 
£(d|0 n ) 



(16) 



We assume uniform prior for the parameters within the 
given ranges. We have made tests with various ranges, for 
the results presented here we have adopted the following: 

n mo e (o,i), n xi e (0,2), a X2 e (0,2), Wl e (-20,5) 

and W2 G (—20,5). We have considered also the special 
case of flat geometry case, for which X1 £ (0, 1 — fi mo ) 
is set, so that £l X2 = 1 — Q mo — Sl Xl > 0, assuming the 
energy density of all the components are non-negative. 

Although the chains will eventually "burn in" even 
given a random position in the parameter-space, in prac- 
tice we have chosen initial points near the maximum like- 
lihood point so that computational time is saved [Hj]. 

The proposal density have an important effect for good 
convergence and mixing of a chain, and the step provided 
by proposal density should not be either too large or 
small that may lead to slow convergence or mixing (39l . 



|43j, |4J, |45j, M, |43, m M, Mk- The method we use to 
determine the proposal density is a adaptive step size 
Gaussian sampler described in Ref. [46j. We use the 
criterion proposed by Gelman and Rubin (1992) to test 
convergence and mixing of the chains, and after that we 
freeze in the the proposal density [46|, [5l[ . 

We generate six chains for each case we study, and 
about seventy thousands points are sampled in each 
chain. After the burn-in process and thinning the chains, 
we merge them into one chain which consists of about 
10000 points used to generate the probability distribu- 
tion of the parameters. As the two dark energy compo- 
nents are completely symmetric for Monte Carlo, when 
plotting the data we have also made "mirrow" points for 
each point on the chain ((Qi,wi) «-> (f^,^). 



D. Model selection criteria 

We shall compare our two component models with one 
dark energy component models. As the number of pa- 
rameters differ, the simple \ 2 statistic is not always ef- 
fective. We therefore consider several model selection 
criteria introduced in the literature [52l IBlil l54ll . 

The Akaike information criterion (AIC) [55| is defined 

as 



AIC 



-2 ln£„ 



2k, 



(17) 



where C max is the maximum likelihood and k is the num- 
ber of parameters in the model [H, [56|, H3] • The second 
term penalizes models with more parameters. However, 
the size of the data set is not considered by the AIC, thus 
when we have very large number of data points, the re- 
duction on x 2 due to additional parameters would also be 
very large, and one may still mis-select the model using 
the AIC criterion 



52] 



The Bayesian information criterion [5a ] can be written 



as 



BIC = -2 ln£ r 



k IniV, 



(18) 



which includes the penalization of the number of data N. 
Nevertheless, the BIC tends to over penalize the num- 
ber of parameters, given large number of data and is an 
approximation of the Bayesian evidence on some assump- 
tions that may not be valid in practice, and it dose not 
take the full advantage of Bayesian technique [52|, [HJ, [56| 
The Bayesian evidence (BE) of a model M takes the 
form 



BE = / C(d\9,M)p(9\M)d9, 



(19) 



which is just the denominator of equation (fTi|) . The BE 
is the average of the likelihood of a model weighted by 
its prior in the parameter space. It automatically in- 
cludes the penalties of the number of parameters and 
data, so it is believed to be more direct, reasonable and 
unambiguous than the Xmin an d I^s in model selection 
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[5J, 15J, 160, |61|, 122, 163, ISJ, M, M ■ The logarithm of BE 
can be used as a guide for comparing models (Jeffreys 
1 961), and we choose the ACDM as the referenced model 
for comparison: Aln(BE) = ln(BE) mof j e / — 1ii(BE)acdm- 
The strength of the evidence for the model is considered 
according to the numerical value of BE: 



A In BE < 1 
1 < A In BE < 2.5 
2.5 < A In BE < 5 
A In BE > 5 



Weak 
Significant 

Strong to very strong 
Decisive 



(20) 



We use the nested sampling algorithm to compute BE 
[13, [fiB], and the error is about 0.10. 



III. RESULTS 




For two component dark energy models, one can define 
an effective equation of state, 



eff(z) 



(21) 



In Fig. [T] we show the distance modulus fi and equation 
of state w e f f for a few examples of two component dark 
energy model, togather with the ACDM model and an 
one-component dark energy model for reference. As can 
be seen from these examples, the phantom divide line 
(w = —1) is easily crossed, which is not unexpected if 
one component has w > — 1 and another has w < — 1. 
The w e f f and [i curves are very smooth, and the the two 
component model can fit the current data pretty well. In- 
deed, from these curves we may take note of the following 
character of the two component model X 1 X 2 CDM1: the 
effective dark energy equation of state approaches that of 
a constant value at high redshift, then it make transition 
to another value at some lower redshift. XCDM2, the one 
component dark energy model with variable equation of 
state with the usual parametrizations, can not realize this 
easily. Inspecting the fit to the distance moduli curve, it 
is not easy to distinguish between these models from each 
other, or from the ACDM model and the one component 
dark energy model, despite their apparant difference in 
the w e ff(z). The error on w derived from current data is 
still too large (the 1 — a error on w for the one component 
model is plotted as yellow dashed lines in Fig. |T|) . 

We now investigate the fit of the two component mod- 
els more quantatively, with MCMC simulation. We have 
performed our fit with several different combinations of 
data sets. We found that the general features of the re- 
sults are all very similar. Here we present only one data 
set of SN so i + GRB in detail. We plot the PDF of fl m0 
in Fig. H The PDF of Q mQ peaks at 0.37 (0.22) for flat 
(generic) geometry. This value is typical for SN la only 
fit, but greater than those based on the combined fit to 
SN la, CMB and LSS data. As we have mentioned in the 
introduction, the purpose of this paper is to study the 
effect of two component dark energy on low redshift SN 



FIG. 1: The w e f / (upper panel) and luminosity distance mod- 
ulus (j, (lower panel) for four two component models, marked 
with red(solid), green (dot), blue (short dash dot) and cyan 
(long dash dot) lines. All these models have flat geometry 
and f2 x i = 0.5, £l X 2 = 0.1, but with different W\,W2- These 
model gave good (but not the best) fit to the data amongst 
the two component models. For comparison, the ACDM 
model with fl m o = 0.27 and the one component model with 
H m o = 0.41, w — —1.5 are plotted. For the latter, we also 
plot the 1 — a error on w in yellow (dashed) line. In the \i 
plot, we show the SN and GRB data points and error bars in 
black and pink respectively. 



Ia fits, we do not include CMB data. There is some dif- 
ference in the best fit (maximal point of the PDF curve) 
value of Cl m o in the one and two component models, but 
the difference is not too large. 

The PDF of one of the dark energy component, fi x i 
is plotted in Fig. [3] For reference, we have also plot- 
ted the PDF of dark energy density in one-component 
dark energy models. Let us consider the flat case first. 
The PDF has a high plateau at fl x i < 0.65 but declines 
sharply above that. As we have f2 m o + O x i < 1 and 
J] m o > 0, this result is not surprising. There are also 
relatively small peaks and two small troughs on the high 
plateau: the peaks at 0.16 and 0.54, and troughs at 0.02 
and 0.37. We shall study the nature of these in the fol- 
lowing paragraphs. However, the peaks and troughs are 
not too large: the amplitude is about 1/10 of the aver- 
age. For the non-flat model, there is a single peak at 0.2, 
and the PDF declines above the peak, but has a small 
low plateau between 0.45 to 0.65. 

In Fig. [4] we plot the historgram for the distribution 
of the equation of state parameter wi. Again, for ref- 
erence purpose we also plotted the distribution for one 
component dark energy models. For both the flat and 
general geometry case, the peak of the distribution is 
near w = — 1. However, for all these distributions there 
are a long tails at w\ <C — 1, for the non-flat case the tail 
is more extended. 
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FIG. 2: The one dimensional PDF of f2 mo . The red solid line FIG. 4: The PDF of equation of state parameter wi, and one 

is for case of flat geometry, the red dot-dashed line for generic component model are also shown, 
geometry. As a comparison, the PDF of one component model 
with static w are shown, and the blue dashed line is for flat 

geometry, the blue dotted line for generic geometry. i — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — 




FIG. 3: The one dimensional PDF of fi-ci. Red solid line: 
flat; red dot-dashed line: generic. For one component model 
with static w, blue dashed line: flat; blue dotted line: generic. 

In Fig. [5]we plot the distribution of f^i vs. W\. We can 
see that the density peak of the distribution is at f2 x i ~ 
0.55, wi ~ — 1. This corresponds well with the peak in 
the PDF of fi x i and w\. On the other hand, we note 
that the distribution also spans to much negative value 
of u>i with small Q X V the addition of a small fraction of 
very negative equation of state component (phantom) is 
allowed by the data. On the other hand, for these deep 
phantom, large O x j is not allowed, leaving a conspicuous 
blank lower right half. The distribution do extends to 
a "pier" at w\ = — 1, Q, X \ > 0.6. These are similar to 



ACDM in the parameter space. There is also some points 
on the corner with w\ > 0. 

In the above we have looked at the distribution of one 
component in our two components model. Now we can 
also examine how the two components are related. The 
scatter plot of £l x iQx2 is shown in Fig. [6l Since fl x i + 
Q x % = 1 — Sl m o, it is very nature to see that there is a 
linear anti-correlation between £l x i and tt X 2- Apparantly, 
the region with 0.4 < f2 x i+f2 K 2 < 0.7 has greater density 
of points, furthermore, in this region, there is a decrease 
of point density near Q x \ — Q X 2- 

We show the scatter plot and associated contours for 
W1W2 in Fig. [71 Near the w\ — W2 diagonal the PDF 
is highest. Models in this region is very similar to one 
component dark energy. As we already know, the one 
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FIG. 6: The point distribution of energy density for the two 
components. 




FIG. 7: The point distribution and contour map of equation 
of state for the two components. The 1<t(68.3%), 2cr(95.5%) 
and 3<r(99.7%) confidence levels are marked by red, green and 
blue solid lines respectively. 



component dark energy model fits reasonably well the 
SN la data, so this result is not surprising. However, we 
note that points far from the diagonal are also allowed, 
indeed the allowed region is extends to very large negative 
w. These models have one dark energy component which 
has —2 < w < 0, and another component whose w can 
be smaller or much smaller than -1. 

Another way to see how the two components are re- 
lated is to plot their differences, as shown in Fig. for 
|f2 x i — 0^,2 1 v s \wi — u>2|. First, we note that the points 
expands to very high value of \u>i — W2[, showing that 
the two components can be very different. Secondly, 



these high difference points are widely distributed from 
\Vl x i — Sl X 2\ — to 0.5. This tells us that it is not the 
case that one single component always dominates. Nev- 
ertheless, we can see that at |fi x i — 0^,2 1 ~ 0.4 the points 
are indeed somewhat more frequent. 




0.4 0.6 
10,,-QJ 

FIG. 8: The difference of the two components. 



IV. MODEL COMPARISONS 

In this section we compare the two component model 
considered in this paper with several different models. 
We calculate their Xmin> as wen as sever al frequently 
used model selection criteria, namely AIC, BIC and BE. 
The following models are considered: ACDM, XCDM1, 
XCDM2, XiXaCDMl, X 1 X 2 CDM2, XiX 2 X 3 CDM (see 
Table 1). XCDM1 is an one component dark energy 
model with static w; XCDM2 is an one component model 
with variable w parameterized by w — Wq + w a z/(l + z); 
X1X2CDMI is the two component model discussed ex- 
tensively in this paper. We also calculate the statistics of 
other two kind of multi-component models: X1X2CDM2 
is two components model with variable w parameterized 
by W!+w la z/{l + z) and w 2 +w 2a z/(l + z); XiX 2 X 3 CDM 
is a three component model with static w. 

With different data set, the best fit result is also differ- 
ent. Here we consider several different data sets: (1) the 
selected supernovae samples and GRB samples, for which 
detailed result is presented in the last section; (2) the 
selected supernovae sample, but without GRB; (3) 182 
supernovae known as the Gold06 sample (4) The SNLS 
sample, which includes 115 supernovae sample published 
in Ref. [2(|, for this sample the parameters a and f3 are 
joined into our cosmological parameter set to estimate 
the statistics; (5) the ESSENCE data set, which consists 
of two parts: 60 ESSENCE samples and 45 nearby SN la 
published in Ref. [2l|. 
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Using the above data sets, we calculate the x 2 and 
other model selection criteria. Our calculation is done for 
flat and non-flat geometric prior seperately. The results 
are given in Table 2 and Table 3. 

First we look at the flat case. It is clear that with addi- 
tional parameters, the % 2 can always be reduced for mod- 
els with more parameters, hence it not very useful in com- 
paring models with different number of parameters [6?| . 
On the other hand, using different model selection cri- 
teria, the conclusion we draw could differ. The BIC 
strongly penalizes models with more parameters, so mod- 
els with less parameters almost always win. For the AIC, 
the results varies with different data set. For the se- 
lected SNIa data set SN se i (which are selected for the 
best data) and SN sc i + GRB data set, the XiX 2 CDMI 
fares as well as XCDMI (the one component dark energy 
model with fixed w), and better than the XCDM2 (the 
one component dark energy model with variable w). In 
the ESSENCE data set, the XiX 2 CDMl fits almost as 
well as XCDMI and better than XCDM2. In the Gold 06 
and SNLS data sets, however, the one component models 
fit better than two component models. 

Perhaps the most accurate model selection criterion 
considered here is BE. Comparing to ACDM, for most 
data set the two component model X1X2CDMI is slightly 
disfavored using BE (except for the ESSENCE data set, 
for which it is slightly favored). However, the evidence of 
this disfavorance is weak. In comparison, the one com- 
ponent dark energy model XCDMI and XCDM2 are sig- 
nificantly disfavored when compared with the ACDM. 
We also found that XiX 2 CDM2 and XiX 2 X 3 CDM have 
comparable BE. 

While we have reason from both theory (inflation) and 
observation (CMB) to suppose that the geometry of ob- 
servable universe is flat, for generality we also consider 
non-flat geometry. In the non-flat case, again with dif- 
ferent data sets the result are somewhat different. Gen- 
erally, however, the data seems to be less favorable for 
two component dark energy models. For example, with 
SN sel + GRB data set, the XiX 2 CDMl is significantly 
disfavored over ACDM, which is similar to XCDMI, 
while the XCDM2 is only slightly disfavored. However, 
XiX 2 CDM2 and XiX 2 X 3 CDM is slightly favored. We 
suppose that in the non-flat model, the curvature term 
behavior is degenerate with a special form of dark energy 
component, which makes the model less favored by the 
model selection criteria. 



V. CONCLUSION 

In this paper we considered phenomelogical multi- 
component dark energy models, in particular a two com- 
ponent model with static equation of state W\ and u> 2 , 
and examined how good can they fit the supernovae and 
GRB data. We have tested our models with several re- 
cently published data sets. For most of the analysis, we 
have used a high quality data set with 182 SN la sample 



and 27 GRB sample selected from resent observations. 
The MCMC technique is used in our cosmological pa- 
rameter fitting, and to compare different models we have 
tried several frequently used model selection criteria. 

The fitting results indicate that the two component 
models can fit the current data set fairly well. The ef- 
fective equation of state of the dark energy can cross 
the phantom divide line, as would have been expect for a 
model with one component w > — 1 and another w < — 1. 

We have compared the multi-component models with 
the ACDM model and one component models using sev- 
eral model selection criteria, including \ 2 , AIC, BIC, and 
BE. The result varies for different choice of model selec- 
tion criteria and data set. The BIC almost always favor 
models with less parameters. For AIC, we found differ- 
ent results for different data sets. Since the Bayesian 
evidence (BE) is presently considered to be better than 
either AIC or BIC, we have mainly focused on the results 
derived using BE. Using BE, we found that for flat ge- 
ometry the multi-component models can fit as well as the 
ACDM model, and for all of the data sets, which is wor- 
thy of note, the multi-component models fit significantly 
(Aln(BE) > 1) better than one component models. 

The model selection technique has been applied to dark 
energy model studies in several recent works. For the 
most recent data, the ACDM model is favored over mod- 
els with simply parametrized variable equation of state 
models such as XCDMI and XCDM2 [68|,[69|. We have 
reached the same conclusion for these models. Because 
the data set and constraining technique used in each case 
is different (e.g., in Ref. (68|, SNIa Gold 06 data set, to- 
gather with CMB shift parameter and SDSS BAO result 
is used), it is difficult to make a full comparison. We note, 
however, that for XCDMI and XCDM2 they obtained 
Aln(BE) = -1.027,-1.118 (with respect to ACDM) re- 
spectively. We have not used the CMB and BAO data in 
our calculation, but for the Gold 06 data set, we found 
A ln(BE) = —1.60, —1.16 respectively, which are of simi- 
lar strength. 

It is premature to say at this stage whether single or 
multiple component dark energy model is favored, as the 
evidence is still weak. Even if in the future we find that 
multi-component model is preferred in a fit, we note that 
this is only a phenomelogical model based on parame- 
trazation. It is conceivable that dark energy of a single 
physical origin could give rise to the appearance of mul- 
tiple components with different equation of state. Never- 
theless, our investigation show that multiple component 
is allowed by current data. When one discusses the prop- 
erties of dark energy as inferred from observations, this 
possibility should not be neglected. 
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TABLE I: H(z) of the models 



Model 


Flat 


Non-flat 


ACDM 


H*{z) = m[Q. m0 {\ + zf + (1 - m0 )] 


H'\z) = H'o[n m o(l + zf + Oao + (1 - o m0 - o A0 )(i + zf] 


XCDM1 


H 2 {z)=H£[n m0 (l + zf 

+ (l-« m0 )(l + 2) 3(1+t " ) ] 


H 2 {z) = H-£{n m0 (i + zf + 0,(1 + 2)^+^ 
+(1-O m0 -0,)(l + «) 2 ] 


XCDM2 


H 2 (z) = H'i[Q m0 (l + zf 

+(1 - o m0 )(i + zf (1+w+w ^e 3w ^ 1/(1+z) -^] 


H 2 (z) = ^[O m0 (l + z) 3 + 0,(1 + z f^ +w+w ^ 

^.[1/(1+,)-!] + (1 _ fimQ _ + ^ 


X1X2CDMI 


H 2 (z) = # 2 [O m0 (l + zf + 0,i(l + zf {1+wi> 

+(i-n m0 -Oxi)(i + «) 3(1+W2) ] 


H 2 (2> = i? 2 [O m0 (l + zf + 0,i(l + 
+O rf (l + zf^ +w ^ + (1 - O m0 - 0,! - 0, 2 )(1 + zf] 


X1X2CDM2 


ti (z) = H [il m0 (l + z) 
+0,i(l + z ) 3 ( 1 + l "i+ 1 "i») e a " , i<.l 1 /(i+^)-i] 

+ (1 - O m0 - 0,0(1 + «) 3 ( 1 +-2+-2a) 
e 3«>2o[l/(l+«)-l]] 


ii (2) = ii ["mo(l + z) 
+0,i(l + 2 ) 3 ( 1 + l "i+ l " 1 "> e 3u, i"[ 1 /( 1 + z )~ 1 l 

+ 0l 2 (l + z) a ( 1 + -2+-2a) e 3»2a[l/(l + Z )-l] 

+(1 - O m0 - 0,i - 0, 2 )(1 + zf] 


X1X2X3CDM 


# 2 ( 2 ) = /f 2 [o mo (i + *) 3 + n-i (1 + ^j 3 ™' 
+o, 2 (i + 2 ) 3(1+ ™ 2) 
+(1 - o m0 - n«i - n„ 2 )(i + ^) 3(1+ ™ 3) ] 


H 2 (z) = tf 2 [O mo (1 + zf + 0,, (1 + zf<- L +™i> 
+0,2(1 + zf ( - 1+w ^ + 0, ;J (1 + z ) 3 ( 1+u, 3) 
+(1 - O m0 - 0,i - 0,2 - 0, 3 )(1 + zf] 



TABLE II: Flat case 



Data 


Type 


ACDM 


XCDM1 


XCDM2 


X1X2CDMI 


X1X2CDM2 


X1X2X3CDM 


Gold06 
(182) 


Xmin 

AIC 
BIC 
Aln(BE) 


158.749 
160.749 
163.953 


156.584 
160.584 
166.992 
-1.60 


156.405 
162.405 
172.017 
-1.16 


156.453 
164.453 
177.269 
-0.68 


155.052 
167.052 
186.404 
+0.04 


156.430 
168.430 
187.654 
+0.15 


SNLS 
(115) 


Xmin 

AIC 
BIC 
Aln(BE) 


110.033 
117.033 
125.268 


110.998 
118.998 
129.978 
-1.26 


111.008 
121.008 
134.733 
-1.34 


110.704 
122.704 
139.173 
-0.36 


110.110 
126.110 
148.070 
+0.66 


109.707 
125.707 
147.666 
+0.64 


ESSENCE 
(105) 


Xmin 

AIC 
BIC 
A ln(BE) 


104.531 
106.531 
109.185 


103.459 
107.459 
112.767 
-1.11 


103.391 
109.391 
117.388 
-0.91 


99.469 
107.469 
118.085 
+0.54 


98.150 
110.150 
126.074 
+0.98 


97.311 
109.311 
125.235 
+0.94 


SN scl 
(182) 


\min 

AIC 
BIC 
Aln(BE) 


162.508 
164.508 
167.713 


161.267 
165.267 
171.675 
-2.23 


161.196 
167.196 
176.808 
-2.70 


157.332 
165.332 
178.148 
-0.50 


153.754 
165.754 
184.978 
+0.52 


155.512 
167.512 
186.737 
+0.40 


SN scl + GRB 
(209) 


Xmin 

AIC 
BIC 
Aln(BE) 


181.341 
183.341 
186.683 


180.029 
184.029 
190.713 
-2.35 


179.988 
185.988 
196.015 
-2.58 


176.405 
184.405 
197.774 
-0.31 


174.618 
186.618 
206.672 
+0.48 


173.466 
185.466 
205.520 
+0.17 



[69] R. Lazkoz and E. Majerotto, JCAP. 0707, 015 (2007) 
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TABLE III: Non-flat case 



Data 


Type 


ACDM 


XCDM1 


XCDM2 


XiX 2 CDMl 


XiX 2 CDM2 


X1X2X3CDM 






156.448 


156.447 


156.206 


156.354 


155.303 


156.434 


Gold06 


AIC 


160.448 


162.447 


164.206 


166.354 


169.303 


170.434 


(182) 


BIC 


166.856 


172.059 


177.022 


182.374 


191.731 


192.863 




A ln(BE) 


- 


-1.45 


-1.15 


-1.66 


-1.07 


-1.61 




'X.min 


111.006 


111.066 


110.495 


110.641 


109.711 


109.544 


SNLS 


AIC 


119.006 


121.066 


122.495 


124.641 


127.711 


127.544 


(115) 


BIC 


129.986 


134.790 


138.965 


143.855 


152.416 


152.248 




Aln(BE) 


- 


-0.80 


-0.07 


-0.78 


+0.46 


-0.14 




Xmin 


103.591 


99.443 


99.393 


98.464 


98.369 


97.303 


ESSENCE 


AIC 


107.591 


105.443 


107.393 


108.464 


112.369 


111.303 


(105) 


BIC 


112.899 


113.405 


118.008 


121.734 


130.947 


129.881 


Aln(BE) 


- 


-1.10 


+0.29 


-1.07 


-0.14 


-1.30 




2 

'Xmin 


161.225 


158.686 


153.607 


157.693 


153.932 


155.108 


SN scl 


AIC 


165.225 


164.686 


161.607 


167.693 


167.932 


169.108 


(182) 


BIC 


171.633 


174.298 


174.423 


183.713 


185.156 


191.536 




Aln(BE) 




-1.10 


+0.90 


-1.34 


-0.24 


-0.40 




\min 


180.397 


178.631 


178.605 


176.608 


174.446 


173.463 


SN sol + GRB 


AIC 


184.397 


184.631 


186.605 


186.608 


188.446 


187.463 


(209) 


BIC 


191.082 


194.658 


199.975 


203.319 


211.843 


210.859 




Aln(BE) 




-1.32 


-0.11 


-1.27 


+0.08 


+0.03 



